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As technological advances allow us to fabricate smaller autonomous self-propelled devices, it is 
clear that at some point directed propulsion could not come from pre-specified deterministic periodic 
deformation of the swimmer's body and we need to develop strategies to extract a net directed 
Q\ , motion from a series of random transitions in the conformation space of the swimmer. We present 

f*"^ ■ a theoretical formulation to describe the "stochastic motor" that drives the motion of low Reynolds 

number swimmers based on this concept, and use it to study the propulsion of a simple low Reynolds 
CN| , number swimmer, namely, the three-sphere swimmer model. When the detailed-balanced is broken 

and the motor is driven out of equilibrium, it can propel the swimmer in the required direction. 
The formulation can be used to study optimal design strategies for molecular-scale low Reynolds 
number swimmers. 
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PACS numbers: 07.10.Cm, 82.39.-k, 87.19.St 



I. INTRODUCTION 



Biological molecular motors [l| are ingenious nano-scale machines that convert chemical energy into directed me- 
chanical work amid strong thermal fluctuations. With the current miniaturization trend in technology, one naturally 
wonders if it is possible to synthesize devices with similar functionalities (2j. In particular, it is desirable as a first 
step to design autonomous small scale swimmers, which could later on be steered by coupling to a guiding network 
Q or system. These swimmers could be used in carrying cargoes or stirring up fluids at small scales. 
O . There is a significant complication in designing swimmers at small scale as they have to undergo non-reciprocal 
deformations to break the time-reversal symmetry and achieve propulsion at low Reynolds number [3[ ■ While it is not 
so difficult to imagine constructing motion cycles with the desired property when we have a large number of degrees 
■ of freedom at hand — like nature does, for example — this will prove nontrivial when we want to design something with 
only a few degrees of freedom and strike a balance between simplicity and functionality, like most human -eng ineered 
C^l ■ devices B. Recently, there has been an increased interest in such designs 0, @, 0, H, H E3, EL EH, [H, 0, Il5l. [Lsl. IrH 
' [l8l . [l^ . |2Ct l2ll . WX . [23] and two interesting examples of such robotic micro-swimmers have been realized experimentally 
using magnetic colloids attached by DNA-linkers 0, HBL Among others, a simple swimmer model based on spheres 
connected by arms that do not interact with the fluid [5[ has been recently used for a number of studies including 
scattering of two swimmers [l?], El, collective hydrodynamic coupling of swimmers [H, HcJ, general feasibility of 
various design properties of swimmers [2l[ , and the effect of large cargos on the performance of swimmers [22| . While 
constructing small swimmers that generate surface distortions is a natural choice, it is also possible to take advantage 
of the general class of phoretic phenomena to achieve locomotion — as they become predominant at small scales — as 
recent experimental [ly, [27], [H, [2^] and theoretical [10, HI], H3 works have demonstrated. 

Here we construct a general statistical mechanical formulation for studying low Reynolds number swimmers that 
undergo conformational changes in a stochastic manner pertinent to systems of molecular scale. We attribute transition 
rates to each deformation move or swimming stroke, and calculate the propulsion velocity as a function of these rates. 
Our formulation provides a general prescription on how to construct the relevant portions of the configurational space 
of swimmers, and how to take advantage of the complexities in this space to maximize the efficiency of the swimmer. 
We apply the formulation to the specific example of the three-sphere swimmer model, which yields interesting results. 

The rest of the paper is organized as follows: Section [IT] describes the general formulation of hydrodynamics of 
low Reynolds number swimmers, and it is followed by Section [IIII that is devoted to the statistical mechanics of the 
conformational changes in swimmers. The formulation is applied to the example of three-sphere swimmer model in 
Section IIV1 which is followed by concluding remarks in Section fVl 
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II. HYDRODYNAMICS OF LOW REYNOLDS NUMBER SWIMMING 

We consider a deformable extended body as a system composed of N point-like solid components described by 
their position vectors r a (t). The deformation is related to internal forces exerted between these solid components, so 
that on component a is exerted a net force f a (i), that is in turn applied oton the fluid at r a (t). In our description 
of point-like objects, the hydrodynamic interactions between the objects relate these force to the velocities of the 
components v Q (/j) = r a (t) via the Oseen tensor Hy(r,r') [33[ (roman indices describe spatial components), namely 

<=E M ^//' (!) 



where M^ 3 = 7iy (r Q ,r' 9 ) and summation over repeated roman indices that define the vector components is under- 
stood. The Oseen tensor is the Green function for the Stokes equation with the appropriate boundary conditions 
and its explicit form depends on the problem we are considering. For example, in the simplest case we can treat the 
solid particles as point-like and use the 1/r-type expressions for the off-diagonal components of the Oseen tensor, 
while putting in 1/(671770) for the diagonal components where a is the radius of the particles and r\ is the viscosity of 
water. If necessary, one could also incorporate finite size corrections and the effect of confining boundaries by using 
the appropriate form of the Green function. 
We can now invert equation ([T]) as 

S3 

where is the resistance (friction) tensor that satisfies M^N^ — S a -ySik. 

For a swimmer that is not subjected to external forces, the local and instantaneous forces in the body are subject 
to the constraint 

X> = 0, (3) 

a 

which yields 

$>#V=°- (4) 

a,(i 

Similarly, if the swimmer is not under the effect of a net external torque, an additional constraint applies 

J2(r a - r CM ) x f a = 0, (5) 

a 

where the center of mass (CM) position is defined as r CM = jj^2 a r a - We note that this condition might not in 
general be satisfied, as in the case of a recent experiment on magnetic doublets (25[. When it does hold, however, it 
will introduce additional constraints on the type of motion and conformations that we can prescribe for the system. 
Finally, in sufficiently symmetric systems the torque-free constraint might automatically be satisfied [2lj . 

We now assume that the relative positioning of the body components are prescribed, in a reference frame that 
moves with the average position and orientation of the body. This reference frame, which we call the "body frame" 
hereon, will be constant during one cycle of the deformation in the body. As a result of the deformation, over the 
period of one cycle the object is expected to be displaced by a small amount due to a non-vanishing translational 
swimming velocity and rotated slightly if there is a non- vanishing rotational velocity as well. The combination of the 
displacement and rotation will determine the new position and orientation of the body frame, which will be used in 
the calculation of the next step of the motion and so on. Therefore, in this picture the motions are grouped into 
separate slow and fast degrees of freedom, in the sense that what is happening over one deformation is cycle (fast 
degrees of freedom) will be averaged to determine a net change in the slow degrees of freedom that determine the 
overall average translation and rotation of the swimmer through the liquid along its trajectory. 

We now assume that the relative positioning of the body components Rf^ = rf — are known in the body frame, 
which means that the relative velocities v" — = Rf^ 3 are also known. [34| These relative positions and relative 
velocities need to be prescribed in a such a way that all the necessary geometrical constraints are satisfied, as for 
example, one cannot expect to have arbitrary distances between a number of points that form a body of a given shape. 
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If the shape of the object and the conformational changes are sufficiently symmetric so that the object swims on 
average in a rectilinear fashion, averaging the velocity of any tagged component a over a complete cycle yields the 
total average translational velocity of the body 

(v«) = V trans , (6) 

as the difference between the velocity of the a component and that of the whole body will be in the form of relative 
deformations that average out to zero. For a more general case the object will have a rotational component superim- 
posed with the translational one, and the average velocity of the tagged body component in the body frame will have 
the following form 

<v Q ) = V trans + O rot x ((r Q - r CM )) , (7) 

where f2 rot is the angular velocity vector of the body about the center of mass. We can extract the translational 
velocity as 



^trans 



= i E «> ■ (8) 



N 

a 

and the rotational component of the velocity as 

^^E^U^-rHM^), (9) 

a 

where 

4* = E Mfa ~ r?»))(tf - r£ M )) - <(r? - rf M )> <(r« - rf»)) , (10) 

a 

is the average moment of inertia tensor for the object. 

We can single out the velocity of one of the components, say a = 1, and describe all of the velocities in terms of 
this and the prescribed relative velocities, namely vf = v} + Rf 1 - Putting this back in equation j4]), we find 

^-^E^f' ( n ) 

where Lij = J2 a p • Note that one can also choose to specify the forces/tensions in the links instead of the 
relative velocities. In this case it will be straightforward to modify the formulation and calculate the velocities. A 
more general framework would encompass prescriptions relating stresses and deformations. 

We can write the relative positioning of the components in the body frame as Rf^(t) = Rq^ + u"^(t), where u"^ 
denote the deformations of the body about the average shape described by Rq^. If we assume that the deformations 
of the body arc relatively small, we can expand equation pip in powers of the deformations and obtain an expression 
for the instantaneous velocity of the tagged (a = 1) component of the body as 

vl(t) = £ 4^uf + £ Bg^ufut + £ C£lf^uf ufur +■-, (12) 

a,/3 a,/3,7,<5 a,/3,7,<5,/i,i/ 

where the coefficients A-j^"' 3 , -B-^ Q/37<5 , Cyj^ , etc. are purely geometrical pre-factors (i.e. involving only the 
characteristic length scales describing the shape of the body) . Averaging over a full cycle, the contribution due to the 
linear terms u*^ and the symmetric combinations ii^ujf + iiJ. S u" 13 — d{u^v? k S )/dt vanish. Therefore, to the leading 
order, we find the average swimming velocity as 

ii /A A a ^ lS \ 

^ trans = E «> = \ E B ^ 5 - e?»f) - E B ^ s ( -4h ) ' (13) 

where A_4"^ 7<5 is the area element enveloped by the trajectory of the system in the (uj^jit^ 5 ) space, and B"^ 6 = 

jj J^fj, ^!jfc Q/37<5 - Note that AA^ lS / At is not a complete time derivative, and its average over a a full cycle does not 
vanish. A similar expression can be written for the angular velocity. The averaging here denotes time averaging if the 
conformation of the system is prescribed. If, however, the system undergoes stochastic conformational changes, the 
averaging needs to be performed over the distribution of the various conformations. The formulation needed to carry 
out this step of the calculation is developed in the next section. 
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FIG. 1: A typical conformation subspace describing the 
shape of the swimmer. Distinct conformational states are 
identified and connected to one another when transitions 
are permissible, making a graph. The swimming velocity 
will be determined by the sum of the currents in each loop 
of the graph (denoted by different colors here), weighted 
by the area of each loop, correspondingly [see equation 

GUI. 



III. KINETICS IN THE CONFORMATION SPACE 



Let us now consider a conformation subspace of the system corresponding to two representative deformations u\ 
and v>2 (see figure [T]). Since we aim to model molecular systems, we should take into account the stochastic nature 
of the conformational changes and not prescribe a deterministic trajectory for the deformation of the system. We 
identify distinct conformational states of the system, denoted as S n , and construct a kinetic description where the 
deformations of the system are described by transitions between these states with given rates, assuming that they 
occur one at a time and do not overlap with each other. We denote the probability of finding the system in S n as P n 
and the rate for transition m->nas k nm . These probabilities are normalized as ^2 n P n = 1. Connecting the states 
that have permissible transitions between them with links, we find a graph that characterizes the conformational 
kinetics of the system in each subspace, as seen in figure [1] To every link, we can attribute a probability current 

J<nm> — kmnPn ^nm^mj (14) 

and at stationary state we can impose the continuity of current at every node, namely 

, J<nm> = 0. (15) 



Solving the system of equations, we can find all probabilities and currents, and in particular the currents J (a) running 
through all the loops in the graph (see figure [T]). We can then write 

(^) = E- 4 ( q ) j («)' ( 16 ) 

a 

where A(a) is the area enclosed by loop a in the conformation subspace. Equation (|16[) shows that the contributions 
from the different loops act together analogously to circuits in parallel, and therefore, it will be the fastest route that 
will determine the effective swimming velocity. In each loop, however, the different legs are connected in series, and 
the slowest kinetic leg will control the contribution to the effective swimming velocity from each loop (see the example 
below) . 



IV. EXAMPLE: THREE-SPHERE SWIMMER MODEL 

We now focus on the specific example of a three-sphere swimmer model [f| . We define the conformation space of the 
swimmer using the two variables (111,112) that describe the longitudinal deformation of the two arms of the swimmer. 
We assume that the two arms can be in the two states corresponding to either Ui = or U{ = Si, and transit from 
one to the other in an almost instantaneous fashion. This means that the configuration space of the swimmer will be 
made of four distinct states as shown in figure [21 defined by different values of the pair (m,U2), namely: state A for 
(61,82), state B for (61, 0), state C for (0,0), and state D for (0,^)- We then assign transition rates to the system, 
corresponding to the average rate of opening and closing of the arms along the cycle 

A^B^C^D^A. (17) 

kAB kBC kcD k D A 



(0,<5,) 




(')',.()) 



(0,0) 



FIG. 2: Conformation space of the three-sphere swim- 
mer model. This minimal model involves only one 
loop. The convention is such that a net swimming to 
the right requires the system to make more cycles in 
the clockwise direction than in the counterclockwise 
direction. 



Note that in this simple example there is only one loop in the conformation space graph of the system (figure [2]) . 

We can now calculate the swimming velocity as a function the transition rates. Using the general formulation 
described in Sections |TT] and IIII1 we find 



V = (18) 



where K 



with l\ and li being the undcformcd lengths of the two arms and a being the 

radius of the spheres [2l|. The probability current J is a function of the transition rates, which can be obtained 
from the following straightforward algebra. At steady state, the current conservation equations can be written as 
J = IzbaPa - ^-abPb = kcsPB - k B cPc = k DC Pc - k C oPD = k A DPo - k DA PA-, which provide us with four 
equations for the current and the four probabilities, which are also normalized as Pa + P B + Pc + Pd = 1- Solving 
the system of linear equations, we find 

j _ kADk D ck CB k B A - kA B k B ckcpkDA 

Eroplacc A by B, C, D^AD^OC^CS + k A BkBckcD + kAB kAD kjjc + kADkAB ksC ) 



Equations (| 18[) and ([19]) give the swimming velocity of the three-sphere swimmer |23( . 

From equation (|19|) it is clear that if detailed balance holds, then J vanishes as the numerator is zero. Using the 
average steady state current, we can deduce the average period of completing one full cycle of the motion along the 
A^B^C^D^A loop as 

T = J-\ (20) 

We can gain a useful insight by looking at the particular limit where the forward rates are all much higher than the 
corresponding backward ones (k B A 3> kAB, etc.). In this limit, we find 

^ = ^ad + k D c + k c x B + k B \, (21) 

which means that the period for a full cycle is the sum of the time intervals needed to complete each leg of the cycle. 

As another example, we can assume that all of the equilibrium fc^'s arc equal to 1 (for simplicity), and that by 
external action only one of them is modified as k B A = 1 + e. In this case, one can show that equation (|19l) yields 

J = , (22) 

16 + 6e V ' 

which leads to a velocity proportional to the perturbation for small values of e and independent of it if the perturbation 
is very large. The linear dependence can be easily understood for a system that is only slightly driven out of 
equilibrium, and the saturation at large perturbations is because the cycling will then be limited by the other three 
unperturbed transitions. In general, one can see that the slowest leg of the reaction controls the average rate of full 
cyclic motion, which suggests the interpretation that in each loop the different legs are connected in series, in analogy 
to circuits. 



V. CONCLUSION 



We have presented a general formulation that can be used in studying the swimming of a small object that undergoes 
stochastic deformations. The program to follow to this end has two stages: (1) treat the deformations as prescribed 
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and follow the hydrodynamic formulation of Section |TT] to calculate the average swimming velocity in terms of the 
relevant deformation variables. (2) Construct the conformation space of the system based on the deformation variables 
and follow the statistical mechanical description of Section IIIII to work out the contributions to the net swimming 
velocity by various modes of swimming defined as loops in the conformation space. We found that a useful circuits 
analogy can be invoked to describe the efficiency of the swimming, with two notable features: (1) the different modes of 
swimming can be effectively considered to act in parallel, which means that their contributions will be independently 
added to each other to yield the net swimming velocity and therefore the fastest route will be the dominant mode of 
swimming controlling the velocity. (2) In each loop, the different kinetic legs could be considered as acting in series 
with respect to one another, which means that the slowest kinetic leg will control the net contribution to the velocity 
by the loop. 

The formulation also allows us to study the effect of an external force or load on the performance of swimmers. 
External forces both add to the hydrodynamic drag and also affect the performance of the swimming strokes as 
activated moves, as the deformations will involve doing work against or being helped by forces endured by the arms. 
These forces will modify the transition rates, and their effects can be readily accommodated by using the force- 
dependent rates in the kinetic formulation. This effect has been studied for the three-sphere swimmer model, which 
has revealed that the performance of the motor strongly depends on where the force is exerted [23[ . This shows that 
for such small swimmers, the concept of a generic force-velocity response breaks down, which might have interesting 
implications for designing molecular swimmers. 
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